Self-Focusing Dynamics of Coupled Optical Beams 
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We theoretically and experimentally investigate the mutual collapse dynamics of two spatially 
separated optical beams in a Kerr medium. Depending on the initial power, beam separation, 
and the relative phase, we observe repulsion or attraction, which in the latter case reveals a sharp 
transition to a single collapsing beam. This transition to fusion of the beams is accompanied by an 
increase in the collapse distance, indicating the effect of the nonlinear coupling on the individual 
collapse dynamics. Our results shed light on the basic nonlinear interaction between self-focused 
beams and provide a mechanism to control the collapse dynamics of such beams. 
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Nonlinear collapse to a singularity and the accompa- 
nying richness of physical effects is observed in several 
fields of physics, such as optics hydrodynamics [5], 
and plasmas [3]. Self focusing and collapse of optical 
beams with a power above a critical power Per [4] in 
nonlinear Kerr media was first predicted in the 1960's 
[si. Over the past decade, with the advent of high peak 
power ultra-short pulsed lasers, self- focusing dynamics 
of optical beams have revealed a remarkable richness of 
spatial and temporal nonlinear phenomena. These in- 
clude observation of the universal self-similar spatial col- 
lapse profile known as the Townes profile [6] , multiple fil- 
amentation , self steepening and pulse splitting [1, Q , 
multiphoton ionization, and supercontinuum generation 
[l,[lo5, along with saturation and plasma generation that 
typically arrest the collapse [Tlj . 

Related optical beam interactions within the context 
of spatial solitons [HI, have drawn considerable in- 
terest in recent years. Spatial solitons in ID Kerr me- 
dia are shown to interact in a particle-like elastic man- 
ner, where the number of solitons and the correspond- 
ing directions and propagation velocities are conserved 
[li^llJ]- Depending on the relative phase, attractive and 
repulsive forces and power transfer are observed between 
interacting solitons [l5|. In saturable nonlinear media, 
which can support (2+l)D solitons, phenomena such as 
soliton fusi on, fission, annihilation, and spiraling occur 
16l . 17 . 3, The interaction between incoherent soli- 
tons, where the medium responds non-instantaneously, 
exhibits similar effects as those observed with coherent 
solitons in nonlinear saturable media [20]. 

In the regime of optical beam collapse in a Kerr 
medium, where self focusing is dominant and beams do 
not maintain their spatial distribution, only a few theo- 
retical studies of beam interactions have been reported 
[21I , [22I . [23I . [24 1, and to the best of our knowledge no 
experiments have been performed. While in these ini- 
tial theoretical studies several qualitative trends, such as 
repulsion, attraction, and fusion of two beams, were iden- 
tified, the detailed dynamics and the transition to fusion 
of two beams, especially when each beam has a power 
near Per, has not been explored. 



In this Letter, we investigate both theoretically and 
experimentally the spatial collapse dynamics of two co- 
herently coupled beams in Kerr media. We observe re- 
pulsion and attraction between the collapsing beams, and 
a sharp transition to fusion of the beams, which is depen- 
dent on their initial power, spatial separation and relative 
phase. We further show that this transition, accompa- 
nied by a peak in the collapse distance, can be exploited 
to control and manipulate the collapse dynamics of two 
coupled beams. Our results shed light on the basic non- 
linear interaction between self-focused collapsing beams 
and are applicable in different scenarios, including that 
of multiple filamentation and collapse of complex multi- 
lobe beams such as necklace beams [25]. 

To investigate the propagation and collapse of two spa- 
tially separated beams in Kerr media, we numerically 
integrate the scalar (2+l)D nonlinear Schrodinger equa- 
tion (NLSE), neglecting dispersion and high-order terms. 
These assumptions are reasonable as long as the disper- 
sion length is longer than the nonlinear and diffraction 
length scales. We express the (2+l)D NLSE in normal- 
ized units as 



Lnl 



(1) 



where is proportional to the amplitude A{x^ z) = 
Ao'0(/i, i^, C) of the envelope of the electric field, ( = 
z/Ldf is the normalized coordinate in the propagation 
direction, /i = x/ro and u = y/ro are the normal- 
ized transverse coordinates, ro is the characteristic ra- 
dius of the input beam, = d'^ /dji^ + jdv^ is 
the transverse Laplacian, L^j = /cro/2 is the diffrac- 
tion length, Lnl = ^o^2c|Aop/A is the nonlinear length, 
k = 27rno/A is the wave number, A is the vacuum wave- 
length, no is the linear index of refraction, and n2 is 
the nonlinear index coefficient. When the total power 
P = {noc/27r)\Ao\'^ J \tlj\'^djj.du in the beam is greater 
than Per = a(A^/47rnon2), where a is a constant that 
depends on the initial spatial profile [4] , the self- focusing 
term dominates and the beam undergoes collapse. 

In the present study, we assume an initial field com- 
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posed of two spatially separated, parallel, Gaussian 
beams of the form V^(/i,z^,0) = C{e-[(^-^o/2)'+^'] ^ 
gi(^g-[(/i+Ao/2)2+^2]^^ where Aq is the normalized initial 
spatial separation of the beams, is the relative phase 
between the beams, and C is a normalization constant 
that depends on the total power. To exclude cases where 
the two beams are initially unresolvable and hence col- 
lapse as a single beam, we assume Aq > 1. Since a col- 
lapsing Gaussian beam is known to form a Townes distri- 
bution carrying one critical power , we limit our inves- 
tigation to the regime of powers 2Pcr < P < 10 Per- This 
corresponds to powers Per < Pg < ^Pcr for each beam, 
thus enabling their individual collapse, but restricts the 
collapse distance to be larger than 0.25Ldf. 

For a fixed total power of 2. 5 Per, various initial separa- 
tions Aq, and for a relative phase of either or tt (in- and 
out-of-phase cases), we integrate Eq. ([T]) with the split- 
step method and obtain the intensity distribution as a 
function of (. Just before collapse, when the peak inten- 
sity reaches 50 times the initial peak intensity, we record 
the C value and the separation between the two collaps- 
ing peaks. The distance between the collapse points Acoi 
and the collapse distance (coi as functions of the initial 
separation are shown in Figs. 1(a) and 1(b), respectively. 
When the separation is large, for both the in- and out- 
of phase cases, the collapse dynamics of each beam are 
essentially independent of the other beam, and the dis- 
tance between the collapse points is equal to the initial 
separation. For the in-phase case, as the initial sepa- 
ration is decreased the collapsing beams are attracted 
towards each other, and their separation decreases with 
respect to the initial separation. At some critical separa- 
tion (Aq ~ 1.9) a sharp transition to one collapsing beam 
is observed. This transition to fusion of the two beams is 
accompanied by more than a 40% increase in the collapse 
distance with respect to independent beam collapse. The 
lengthening of collapse distance represents a critical slow- 
ing down, associated with the phase transition between 
independent collapse and fusion. For the out-of-phase 
case the behavior is drastically different. When the initial 
separation is decreased the collapsing beams repel each 
other and do not undergo fusion to a single collapsing 
point. This resembles the well-known repulsion observed 
between two anti-phased spatial solitons [12|]. Here the 
collapse distance is decreased with regard to independent 
beam collapse. 

We further investigate for the in-phase case the depen- 
dence of the above transition on the total power. Fig- 
ure 1(c) shows the distance between collapse points as 
a function of initial separation for different total pow- 
ers. As the total power is decreased, approaching 2Pcr 
(i.e. one critical power per beam), the critical separation 
at which the transition occurs increases, and the transi- 
tion becomes steeper. With a total power of 2.02Pcr, the 
critical separation is almost three times that for a total 
power of lOPcr- As P — > Per, the distance to collapse is 
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FIG. 1: (Color online) Numerical calculations of the collapse 
dynamics of two coupled Gaussian beams in Kerr media, (a) 
Spatial separation between collapse points and (b) collapse 
distance as a function of initial separation between the beams, 
for and tt relative phase (total power 2.5Pcr); (c) same as 
the in-phase case in (a) for various input powers. 



extended, allowing for longer lengths of mutual attrac- 
tion and eventually for fusion into one beam (e.g. with a 
total power of 2.5Pcr the independent collapse distance 
is lALdf while with 2.02Pcr it is TL^f, and is more than 
doubled at the transition region to ISLdf). At relatively 
large total powers, the transition occurs at smaller initial 
separations, such that the beams appear to be decoupled 
for most initial separations. Unlike the case of soliton in- 
teractions, these dynamics are governed simultaneously 
by two competing processes: individual beam collapse 
and fusion of the beams. 

To experimentally investigate the collapse dynamics of 
this system, we use the arrangement schematically shown 
in Fig. 2. A pulsed, Ti: Sapphire amplified laser beam, 
with 1-kHz, 90-fs pulses centered at 800 nm, is passed 
through a spatial filter to obtain a nearly Gaussian pro- 
file. This beam is sent into a two- arm interferometer. 
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composed of a 50:50 beam splitter and two retro reflec- 
tors, such that the two spatiahy separated beams can 
be recombined at the same beam sphtter. The separa- 
tion of the two beams is changed by moving one of the 
retro reflectors in the transverse direction. The tempo- 
ral alignment needed for coherently combining two short 
pulses and for controlling their relative phase is achieved 
by changing the length of one arm with a high-precision 
motorized stage. The combined beams emerging from 
the interferometric setup are telescoped down to obtain 
a 0.29-mm beam waist (cjq) and then sent into a 28-cm 
rod of BK7 glass (equivalent to 0.56Ldf). The intensity 
distribution at the input and output faces of the BK7 rod 
are imaged onto a 12-bit CCD camera. In order to mon- 
itor the relative phase of the two combined beams, the 
far-fleld intensity fringes of the input fleld distribution 
are imaged onto an additional CCD camera. 

Initially, one arm is blocked, and the laser power is ad- 
justed so that collapse of a single beam is obtained at the 
output face of the BK7 rod. The energy per pulse is mea- 
sured to be 3.6 /i J, equivalent to 40 MW, which is about a 
factor of 10 higher than expected from the CW numerical 
calculations for collapse after O.bGLdf in BK7 glass. We 
attribute this discrepancy in powers mainly to dispersive 
effects in the long glass rod and flnd that this translates 
into merely a scaling factor for the power and does not 
affect signiflcantly the spatial dynamics under investiga- 
tion. This power, at which collapse of a single beam 
is observed at the output, is held constant throughout 
the experiments. In practice, due to coherent overlap of 
the beams on the beam splitter, the total power slightly 
increased as their separation is decreased. For beam sep- 
arations greater than l.buoo this increase is calculated to 
be less than 32%. 

After setting the power, we unblock the two arms and 
for each initial separation of the two beams we record the 
input and output intensity distributions. This is done for 
the case in which the beams are temporally misaligned 
(i.e. the independent beam case) and for both the in- and 
out-of-phase cases, where the relative phase is zero and 
TT, respectively. Several representative intensity distribu- 
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FIG. 2: Ultrafast experimental arrangement for self focusing 
and collapse of two spatially separated Gaussian beams in 
BK7 glass. 



tions for different initial separations and relative phase 
are shown in Fig. 3. For the in-phase case, the observed 
behavior is in agreement with our simulations in which 
two collapsing points occur at the output when the ini- 
tial separation is large, whereas only one collapse point is 
seen when the initial separation is small. This is a clear 
indication of predicted attraction and fusion between two 
collapsing in-phase beams. For the out-of-phase case, two 
collapsing points are always observed at the output, in- 
dicating repulsion between the collapsing beams, which 
is also in agreement with our theoretical results. 

From the recorded intensity distributions, we measure 
the initial separation between the beams at the input 
face and the corresponding separation between the col- 
lapse points at the output face, and the results are plotted 
in Fig. 4. As expected from the numerical calculations, 
attraction between the beams and a sharp transition to 
a single collapse point is observed for the in-phase case. 
This occurs at a separation distance of 1.75a;o, corre- 
sponding to the numerical calculations with Ptotai ^ 4Pcr 
[see Fig. 1(c)], which is in good agreement with the col- 
lapse distance. For the out-of-phase case, the beams re- 
pel each other, and the separation between the collapse 
points does not decrease below 0.52 mm, which is equal 
to l.ScjQ. 

In order to experimentally test our ability to dynami- 
cally control the collapse of two coupled beams via their 
relative phase, we slowly sweep the time overlap of the 
two beams (near the optimal overlap point), thus chang- 
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FIG. 3: (Color online) Experimental input and output inten- 
sity distributions of two collapsing beams with different initial 
separations for the in- and out-of-phase cases. The Gaussian 
beam waist is 0.29 mm, the total energy for the two beams is 
7.2 /iJ, and the propagation distance in the BK7 glass sample 
is 28 cm. 
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FIG. 4: (Color online) Experimental measurements of the 
spatial separation between collapse points as a function of 
initial spatial separation between two beams for and tt rel- 
ative phase and for the temporally misaligned case. Here the 
initial beam waist is 0.29 mm, the total energy in the two 
beams is 7.2 /iJ, and the collapsing beams are measured after 
28 cm of BK7 glass, which corresponds to the distance where 
independent collapse occurs. 



ing their relative phase periodically. This is performed for 
two different initial separations, the first at a separation 
of 0.41 mm corresponding to IAcoq where the in-phase 
case results in fusion into one collapse point, and the sec- 
ond at a separation of 0.54 mm corresponding to I.SGcjq 
where the distance between the two collapse points is 
significantly affected by the relative phase. In the first 
case the output distribution oscillates between two dis- 
tinct states with either one or two collapse points. In 
the second case, the output distribution always has two 
collapse points, for which the separation increases and 
decreases periodically. When the relative phase is in the 
intermediate state (neither or tt) , the power appears to 
redistribute between the collapse points in an analogous 
way to that observed in ID soliton experiments [13.]. 

In conclusion, we have investigated the collapse dy- 
namics of two spatially separated coupled beams in a 
Kerr medium. For the in-phase case, we experimentally 
observe a sharp transition to fusion of the two collaps- 
ing beams, which depends on the initial separation and 
power. We predict that this transition is accompanied 
by a longer collapse distance with respect to indepen- 
dent collapse. For the out-of-phase case, the beams repel 
each other and collapse independently. Our results of 
two-beam interaction serve as a basis for analyzing more 
complicated scenarios where multiple beams or compos- 
ite multi-lobe beams undergo self focusing and collapse. 
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